Conductance of a quantum wire in the Wigner crystal regime 
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We study the effect of Coulomb interactions on the conductance of a single-mode quantum wire 
connecting two bulk leads. When the density of electrons in the wire is very low, they arrange in a 
finite-length Wigner crystal. In this regime the electron spins form an antiferromagnetic Heisenberg 
chain with exponentially small coupling J. An electric current in the wire perturbs the spin chain 
and gives rise to a temperature-dependent contribution of the spin subsystem to the resistance. At 
low temperature T < J this effect is small, and the conductance of the wire remains close to 2e 2 /h. 
At T S> J the spin effect reduces the conductance to e 2 /h. 

PACS numbers: 73.63.Nm, 73.21. Hb, 75.10.Pq 



Experiments with short one-dimensional (ID) conduc- 
tors (quantum wires) have demonstrated [jj that their 
conductance is quantized in units of 2e 2 /h. The univer- 
sality of this result is readily understood in the model 
of non-interacting electrons. In this approach the con- 
ductance is given by G = e 2 WFi where vf is the Fermi 
velocity of electrons and v — 2/hvp is their density of 
states in one dimension. It is remarkable, however, that 
in most experiments observing the quantization of con- 
ductance the Coulomb interactions between electrons are 
not weak, (nag) -1 > 1. (Here n is the electron density, 
as = eh 2 /me 2 is the Bohr radius, e is the dielectric con- 
stant, and m is the electron effective mass.) Electron- 
electron interaction in a ID system is expected to lead to 
the formation of a Luttinger liquid with properties very 
different from those of the non-interacting Fermi gas. 

The conductance of an infinite Luttinger liquid was 
studied by Kane and Fisher |2(], who found that it does 
depend on the interactions. In particular, in the case of 
repulsive interactions the conductance is below the uni- 
versal value of 2e 2 /h. The discrepancy between the the- 
ory and experiments []J is usually attributed to the 
fact that in order to measure conductance of a quantum 
wire, it has to be connected to two-dimensional Fermi- 
liquid leads. As electrons move from the wire into the 
leads, the interactions between them are gradually re- 
duced to zero. At low frequency to -C Vp/L the conduc- 
tance of such a system is dominated by the leads, and 
the universal value 2e 2 /h is restored 3]. 

The key assumption leading to this result is the appli- 
cability of the Luttinger liquid description of the inter- 
acting electron system in a quantum wire. In this paper 
we show that if the interactions are strong, {nas)^ 1 ^> 1, 
the Luttinger liquid picture remains valid only at expo- 
nentially low temperatures, and study the corrections to 
the quantized conductance at higher temperatures. 

The low-energy properties of an interacting ID electron 
system are most conveniently described by the bosoniza- 
tion approach. In the case of weak interactions the 
Hamiltonian of the system can be presented Q as 



where 



H P = 
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H — H p + H a , 
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Here the fields cf> p , H p and <f> a , H a describe the excitations 
of the charge and spin modes, respectively, and satisfy 
the bosonic commutation relations: [<p u (x),H u i(x')] — 
i8 llv i5{x — x'). Parameters K p , K ai and gu_ are deter- 
mined by the interactions between electrons, u p and u a 
are the velocities of propagation of spin and charge exci- 
tations, and a is a short-distance cutoff. 

An important feature of the Hamiltonian Q is the 
separation of the charge and spin variables 11]. This 
property is preserved even when bias is applied to the 
wire, because the electric potential couples only to the 
charge density, and does not excite the spin modes. As a 
result, the relatively complicated form JSJl of the Hamil- 
tonian H s is not expected to affect the conductance of 
the quantum wire. 

The derivation Q of the bosonized Hamiltonian as- 
sumes weak interactions between electrons. On the other 
hand, the form (Q-© of the Hamiltonian is universal, 
i.e., with the appropriate choice of the parameters K Pt<7 , 
u p ^ ai and gi±^ it is expected to describe the low-energy 
properties of ID electron liquids with arbitrarily strong 
interactions. It will be instructive to obtain the Hamil- 
tonian CQ-Q from the model of fermions with strong 
Coulomb interactions, and to estimate the values of the 
parameters entering Eqs. (0), J3Jl. 

At low density n <C a^ 1 the potential energy of the 
electron repulsion ~ e 2 n/e is much larger than their 
kinetic energy ~ h 2 n 2 /m. Thus to first approxima- 
tion one can view the electrons in a quantum wire as 
a Wigner crystal of particles repelling each other with 
strong Coulomb forces. Assuming that the Coulomb in- 
teraction is screened at large distances by a metal gate 
parallel to the wire, the density excitations of the Wigner 
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crystal (plasmons) are acoustic waves. Thus the low- 
energy density excitations of this system must be de- 
scribed by the Hamiltonian of the form The plas- 
mon velocity in a one-dimensional Wigner crystal at a 
distance d from the gate is u p — U e 2 n/mC, where 
C = e/[2ln((nd)] is the capacitance per unit length be- 
tween the crystal and the gate [f| . (Here C is a numerical 
coefficient determined by the geometry of the gate [lij.) 
Parameter K p can be found by comparing the second 
term in the integrand of Eq. J2J with the charging energy 
per unit length Eq — e 2 Sn 2 /2C, where 5n(x) is the devi- 
ation of the density of electrons from its mean value. The 
bosonization procedure used in the Hamiltonian J3J iden- 
tifies Sn = —^-d x (j) p . Comparing Ec with the second 
term in Eq. one then finds K p — (nh/2) \JnC / me 2 . 
Substituting the above expression for C, we summarize 
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where the parameter vp = nhn/2m is defined as the 
Fermi velocity in a gas noninteracting electrons at den- 
sity n. As expected, at (nas)^ 1 3> 1 strong Coulomb 
interactions result in K p <C 1 and u p 3> vf- 

In the limit of strong coupling the energy of a Wigner 
crystal does not depend on the spins of electrons. In- 
deed, electrons localized near their equilibrium positions 
can be viewed as distinguishable particles, and the anti- 
symmetrization of the wavefunctions does not affect the 
energies of the eigenstates. This picture is violated if one 
allows for the possibility of overlap of the wavefunctions 
of neighboring electrons. The overlap is due to the pos- 
sibility of tunneling through the potential barrier e 2 jer 
separating the electrons. To first order in tunneling, the 
overlap occurs only between the nearest neighbors, and 
one expects the coupling of the spins to be described by 
the Hamiltonian of a Heisenberg spin chain: 
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Since the ground state of a system of interacting fcrmions 
in one dimension cannot be spin-polarized |fjj], the ex- 
change must have antiferromagnetic sign, J > 0. 

A careful evaluation of the exchange constant J 
presents a challenging problem of many-body physics, 
which is beyond the scope of this paper. A crude esti- 
mate of J can be obtained by calculating the amplitude of 
tunneling through the Coulomb barrier e 2 / ' er separating 
two neighboring electrons in the WKB approximation. 
Placing the turning points at r — ±n _1 and using the 
reduced mass m/2, one finds 



J ~ Ep exp I — 



V 



(6) 



with the numerical coefficient r\ = n. (Here Ep = 
Tr 2 h 2 n 2 /8m is defined as the Fermi energy of a non- 
interacting electron gas of density n.) A more careful 



calculation will likely result in a different value of r/, but 
the exponential dependence of J on nas will remain. 

The Hamiltonian of the Heisenberg spin chain © can 
be rewritten in terms of spinless fermion operators ai and 
a\ with the help of the Jordan- Wigner transformation 



Sf = ajai 



~ Sf 



iSf = a\ exp 



I-i 

nrJ2 a ] a 3 | ■ ( 7 ) 



To study the low-energy properties of the spin chain, one 
can bosonize the fermion operators ai and a[. As a result 
the Hamiltonian H a takes the form ©, see, e.g., Ref. 0. 
The velocity u a of the spin excitations is easily deduced 
from the spectrum of the isotropic Heisenberg model, 



ttJ 
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The sin-Gordon perturbation in Eq. ((3J) is marginally ir- 
relevant, i.e., the coupling constant gi± scales to zero at 
low energies. At the same time the parameter K a scales 
to 1 , as required by the SU(2) symmetry of the problem 

i, m 

It is important to stress that the bosonized form J3J of 
H a is only appropriate at low temperatures, T <C J ■ In 
the following we will also be interested in the tempera- 
ture dependence of the conductance at T ~ J, and thus 
we will use the form (JSJ). Interestingly, the dynamics of 
spin and charge modes are still completely separated, as 
the operators J2J and (J5J commute. We now show that 
this symmetry is violated when the wire is connected to 
Fermi-liquid leads. 

Following Ref. 0, we model the leads attached to 
the quantum wire by two semi-infinite sections of non- 
interacting electron gas. To this end we assume that the 
ID electron density n(x) = n near x = 0, and gradually 
grows to a very high value a^ 1 at x — > ±oo. 

Assuming that the length scale L of the dependence 
n{x) is large compared to the distance between electrons, 
one can neglect the backscattering caused by the inho- 
mogeneity, and describe the charge excitations by the 
bosonized Hamiltonian J2J) with position-dependent pa- 
rameters u p {x) and K p (x). In addition, the coupling con- 
stant J in the Hamiltonian JSJ now depends on I due to 
the density dependence © of the exchange interaction. 
The exchange constant J [I] is determined by the density 
n(xi) at the position of the Z-th electron. 

Now let us consider the quantum wire in the regime 
when an electric current 1 = 1$ cos uit passes through 
it. We will be interested in the dc limit u> — * 0, and 
can thus assume that all electrons move in phase. The 
charge transferred through any point in the wire is Q = 
Iqw^ 1 smuit, so at moment t the Z-th electron has shifted 
to the position Z + io(ew) -1 sina>£. Thus the Hamiltonian 
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of the spin chain takes the form 

H r7 =Y / J[l + q(t)]S l -S l+1 , q(t) = —mjxut. (9) 

l 

One can view current I = eq as an excitation of the 
charge mode <p p and substitute the appropriate relation 



?(*) = — 0,(0, t) 

7T 



(10) 



into the Hamiltonian @. Therefore the spin modes are 
coupled to the charge ones and should affect the conduc- 
tance of the wire. 

In the following it will be more convenient to treat 
the current I(t) — Iq cos Lot as an external parameter. 
This approach corresponds to the experiment with the 
wire connected to a current source. We will evaluate the 
energy W dissipated in the device in unit time in the limit 
of small Iq and lo. The dc resistance will be found from 
the Joule heat law W = \IqR- The advantage of this 
approach is that the spin and charge modes are coupled 
only through the external parameter I(t), so W will be a 
sum of two independent contributions of the charge and 
spin modes. Thus the resistance of the wire is the sum 
of two terms R = R p + which can be evaluated by 
considering the Hamiltonians H p and H a separately. 

The result G = 2e 2 /h for the conductance of a quan- 
tum wire found in Ref.|3j amounts to the calculation of the 
resistance R p , as the spin modes were ignored. Thus one 
expects to find R p — h/2e 2 . Let us outline the deriva- 
tion of this result in the approach where the current I(t) 
through the wire is fixed. Then the charge subsystem 
is described by the Hamiltonian (J2J with the boundary 
condition ^10(1. It will be convenient to transform the 
variables <f> p (x) — > 4>p( x ) + 7r <?(t)/V / 2 J i- e -i to apply to the 
Hamiltonian a unitary transformation 

U = cxp (-i^^- J ILp(x)dx \ ■ (11) 

Upon this transformation the boundary condition be- 
comes time- independent, <p p (Q) = 0, and the Hamilto- 
nian transforms into H p = U^H p U — ihU^dtU. The sec- 
ond term is a time-dependent perturbation that excites 
plasmons of very low frequency lo. The wavelength of 
these plasmons is large, u p /lo 3> £, and thus one can 
replace u p and K p by their values at x — > oo; in par- 
ticular, K p (po) — 1, as riooOB ^> 1. The transformed 
Hamiltonian H p is conveniently presented in terms of the 
operators 



l>i. = / W( /■•.(•) >in A\r | — — -4> P {x) H -j=Y\ p {x) j dx. 



Then H p takes the form 



destroying plasmons with frequency lo^ — u p (po)\k\. 

(12) 



h p = h 



LJkblbk+i^- =g dk. 



At low frequency lo <C T the time-dependent perturba- 
tion leads to both emission and absorption of plasmons 
with k — ±lo/u p (oo). The energy dissipated into plas- 
mon excitations in unit time is easily found by means 
of the Fermi golden rule, resulting in W = \Iq (nfife 2 ). 
Thus we conclude that R p = h/2e 2 . 

To find i?o- one has to perform a similar calculation 
with a more complicated Hamiltonian Performing 
the Jordan- Wigner transformation 0, we rewrite it as 



(13) 



In the absence of the external magnetic field the average 
spin per site (Sf ) = 0. In terms of the Hamiltonian 113|) 
it means that the fermionic band is half- filled, (ajai) = h. 
The exchange J [I] takes its lowest value J at the center 
of the wire and grows up to J [00] ~ Ep in the leads. 

The Hamiltonian i|13|) can be easily treated if one ne- 
glects the interaction term in its second line. (This cor- 
responds to the XY model of a spin chain, in which the 
^-component of exchange vanishes.) Then the Hamilto- 
nian describes a tight-binding model of non-interacting 
fermions with the bandwidth 2 J[l] varying gradually be- 
tween the small value 2 J in the wire and a large value ~ 
Ep in the leads. The fermions with energies |e| < J pass 
through the constriction without backscattering, whereas 
the ones with |e| > J are reflected. The time depen- 
dence J [I + q(t)] can be interpreted as a slow movement 
of the constriction with respect to the Fermi gas, and 
leads to the dissipation of energy. For non-interacting 
fermions the calculation of W and the respective resis- 
tance is rather straightforward, and we find 



R? Y = 



1 



2e 2 e J / T + 1 ' 



(14) 
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At T -C J this result is exponentially small, R^' oc 
e~ J / T , because at low temperature most of the quasipar- 
ticles pass the constriction without scattering. Only an 
exponentially small fraction of the excitations are scat- 
tered at the constriction and acquire the energy from 
the driving field. The resistance saturates at T/J — > 00, 
when all the quasiparticles are backscattered. 

The result (|14|l remains qualitatively correct for the 
full model (|13|l . In particular, at J -C T the transport of 
spin- fermions through the constriction is still suppressed, 
and R a saturates. In order to find R a at T/J — > 00, one 
can notice that the bosonized model 10 is still applicable 
in the leads, where J ~ Ep ^> T. Then the absence of 
spin transfer through the constriction can be expressed 
as a hard- wall boundary condition upon <f> a . 

At q(t) = this boundary condition may be presented 
without loss of generality as <j) a {Q) = 0. At non-zero q(t) 
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the same condition <ft a = is imposed at point I = —q(t). 
The half-filling condition for the Hamiltonian (|13fl means 
that q(t)/2 fermions have passed through the wire at mo- 
ment t. In analogy with the calculation of R p , the bound- 
ary condition can be interpreted as application of a fixed 
current q(t)/2 of the spin-fermions. The bosonization 
procedure leading from Hamiltonian i|13|) to Eq. J2J) ex- 
presses the current of spin-fermions as d t 4> a /^V2. Thus 
the appropriate boundary condition for the bosonized 
spin Hamiltonian J3J is ^^ CT (0, t) — q(t). Note, that 
this boundary condition is equivalent to Eq. (|10|l . and 
the Hamiltonians (0 and © coincide in the leads, where 
K p = K a = 1 and g\± = 0. Therefore one can complete 
the evaluation of R a by repeating the above calculation 
of R p , and we conclude that R a = h/2e 2 . The conduc- 
tance of the wire (R p + reduces to e 2 /h. 

It is worth mentioning that the same boundary condi- 
tion for the Hamiltonian J3J appears even at T < J if a 
sufficiently large magnetic field B is applied. Indeed, if 
the Zeeman splitting is large compared to both T and the 
coupling J, the electrons in the wire are completely spin 
polarized, (Sf) = — i. In terms of the Hamiltonian 
this is interpreted as lowering of the chemical potential 
below the bottom of the fermionic band in the wire, so 
that (ajai) = 0. Thus the two leads are now separated 
by a barrier created by the narrowing fermionic band in 
the wire. The barrier is centered at I = —q(t), and, upon 
bosonization, imposes the same boundary condition upon 
4> a as in the case of B = and T > J. We therefore con- 
clude that in a polarizing magnetic field R a = h/2e 2 , 
and, as expected, the total conductance G — e 2 /h. 

At low temperature T <C J the model (|13[) to first 
approximation can be bosonized to the form @ with 
position-dependent parameters u a , K a and gi±. At T — > 
we have K a = 1 and gi± = 0, i.e., the Hamiltonian © 
can be viewed as the bosonized version of the XY model. 
Then we conclude from Eq. I|14|) that at zero temperature 
R a = 0, and the conductance of the wire remains 2e 2 /h. 

At finite length of the wire L the parameters of the 
Hamiltonian @ do not reach their limiting values K a = 
1, gi± = 0. However, the resulting corrections to R a re- 
main small in 1/nL even at T ~ J. A more interesting 
correction appears due to the fact that the bandwidth 
of the Hamiltonian (|13fl is finite, which is not accounted 
for accurately by the bosonization procedure leading to 
Eq. 10. If the wire is long, nL ^> 1, the spin chain © 
is nearly uniform at each point Its low energy exci- 
tations are conveniently classified Q in terms of spinon 
quasiparticles with spectrum e — \ J [I] ship, where p is 
the wavevector in the lattice model. At low temperature 
T <C J typical excitations have energies e ~ T small com- 
pared to the bandwidth in the center of the constriction. 
These spinons pass through the wire without scattering 
and do not interact with the driving field. On the other 
hand, a small fraction of the spinons have energies ex- 
ceeding ttJ/2. These excitations are not supported by 



the spin chain at the center of the wire, and thus are 
reflected by the constriction. Similar to the case of the 
XY model, the reflected particles acquire energy from the 
driving field and contribute to the dissipation. The den- 
sity of such spinons is exponentially small, and we get 



R a ~ i? exp y-yfj ■ (15) 

This treatment does not enable one to evaluate Rq. The 
analogy with the XY model suggests i?o ~ h/e 2 . 

Experimentally, corrections to the quantized conduc- 
tance showing activated temperature dependence consis- 
tent with H15J1 have been reported in short wires jjj. At 
higher temperatures the conductance tends to saturate 
at 0.7 x (2e 2 /h) instead of e 2 /h. Quantization of con- 
ductance at e 2 /h in the absence of magnetic field was 
reported in longer wires at low electron density [|J. 

To summarize, in the regime of strong interactions, 
nas <C 1 , the electrons in a quantum wire form a Wigner 
crystal with exponentially small antifcrromagnetic spin 
exchange J. At T <C J the conductance remains 2e 2 /h 
up to exponentially small corrections, Eq. (|15|) . At 
J <S^T the conductance drops to e 2 /h. Remarkably, this 
result does not assume spontaneous spin polarization in 
the wire 
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